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It is a puzzle why the densities of dark matter and dark energy are nearly equal today when they 
scale so differently during the expansion of the universe. This conundrum may be solved if there 
is a coupling between the two dark sectors. In this paper we assume that dark matter is made of 
cold relics with masses depending exponentially on the scalar field associated to dark energy. Since 
the dynamics of the system is dominated by an attractor solution, the dark matter particle mass is 
forced to change with time as to ensure that the ratio between the energy densities of dark matter 
and dark energy become a constant at late times and one readily realizes that the present-day dark 
matter abundance is not very sensitive to its value when dark matter particles decouple from the 
thermal bath. We show that the dependence of the present abundance of cold dark matter on the 
parameters of the model differs drastically from the familiar results where no connection between 
dark energy and dark matter is present. In particular, we analyze the case in which the cold dark 
matter particle is the lightest supersymmetric particle. 
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1. Introduction. Combined analysis of cosmological 
observations such as cosmic microwave radiation (CMB) 
anisotropics [1], cluster baryon fraction [2] and super- 
novae la [3], give increasing support for the so called 
'cosmic concordance' model, in which the universe is flat 
(f2tot = l)j an d made for one third of non-relativistic 
dark matter (DM) , and for two thirds of a smooth compo- 
nent, called Dark Energy (DE). DM is commonly associ- 
ated to weakly interacting particles (WIMPs), and can be 
described as a fluid with vanishing pressure. DE, being 
responsible for the accelerated expansion of the universe, 
has to be identified with some more exotic component 
with negative pressure, such as the cosmological constant 
or a scalar field with a proper potential [4]. A common 
assumption is that there is no interaction between DM 
and DE. This implies that the ratio between their energy 



densities evolves in time as Pdm/pde 



-,3W 



a being 



the scale factor of the universe and W the DE equation 
of state (W = —1 for a cosmological constant, W < —1/3 
to get acceleration). As a consequence, the cosmological 
constant or the scale of the effective potential have to be 
carefully tuned in order to get the observed proportions 
of DM and DE right today. 

A possible solution to this 'cosmic coincidence' prob- 
lem may be obtained by introducing a coupling between 
DM and DE. More specifically, if DE is identified with 
a scalar field cp with an exponential potential and the 
DM particle mass also depends exponentially on <p, the 
late time behavior of the cosmological equations gives 
accelerated expansion and, at the same time, a constant 
DM/DE ratio [5]. This behavior relies on the existence 
of an attractor solution, which makes the late time cos- 
mology insensitive to the initial conditions for DE and 
DM. The possibility of varying mass particles (VAMPs) 
as DM candidates was considered previously in [6], but 



with different functional dependences not leading to a 
solution of the cosmic coincidence problem. 

In this letter, we will explore the impact of this sce- 
nario on the usual cold DM picture, which is based on 
the freeze-out at an early epoch of some non-relativistic 
species (see, for instance, [7]). In the standard scenario 
with constant mass particles the late time DM energy 
density is just the one at the freeze-out rescaled by the 
volume expansion. Instead, in the present VAMP sce- 
nario the attractor prompts a mass change in order to 
have the right DM/DE ratio today. As we will see, be- 
sides solving the cosmic coincidence problem, this will 
significantly enlarge the parameter space allowed by the 
observed abundances. 

2. VAMPs and the coincidence problem We will con- 
sider as DM candidate a particle \ of mass M, depending 
exponentially on the DE field ip, 

%)=Ie^, (1) 

where <p is expressed in units of the Planck Mass M p . 
The scalar field has an exponential potential 

V(p) = W i (2 ) 

with A, /3 > 0. If the DM particle is non-relativistic, its 
energy density is also (^-dependent, since it is given by 



(3) 



where Y x = n > 



/n-y is the number density of DM particles 



relative to the density of photons, n 1 
o 
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n~ =411 cm 



with 

3 , and a/ao the scale factor of the universe 
relative to its value today. 

As a consequence, the scalar field is subject to an ef- 
fective potential given by the sum of Eqs. (2) and (3). 
The equation of motion is then 
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\ \Px + Pb + 2V+-p Iad J<p' = -pV + \p x , (4) 

where primes denote derivatives with respect to £ = 
log(a/ao), Pb and p ra d are the energy densities in baryons 
and radiation, respectively, and we have assumed a spa- 
tially flat universe. Since we are interested in the late- 
time behavior, we can assume pb, Prad *C Pxi P>p- 111 this 
limit it is easy to sec that there is a solution such that 
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, (5) 



which is an attractor in field space if (3 > (— A + 
VA 2 + 12) /2. From the moment the attractor is reached, 
the energy densities in DM and in DE evolve at a con- 
stant ratio depending only on A and j3, thus solving the 
cosmic coincidence problem. 

A closer inspection at the solution (5) using Friedman's 
equations reveals that 
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with W 
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(6) 



that is W is negative and may lead, if W < —1/3, to 
an accelerated expansion of the universe. To understand 
how it is possible to get both acceleration and constant 
ratio between DM and DE one may look at the scaling 
behavior of the energy densities on the attractor (5), 
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(7) 



The tp dependence of the DM mass modifies the usual 
scaling e -3 ^ of non relativistic matter. Since the mass 
increases with the expansion, it corresponds to an effec- 
tively negative W, even though DM is still a pressureless 
fluid made up by non-interacting particles. 
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FIG. 1. Time evolution of the relative abundance of 
different species, expressed as fractions of the critical 
density. The corresponding parameter configuration is: 
/3 = 1.88, A = 0.82, V = 10 2 GeV 4 , M s = 500 GeV, 
Mf = 105 GeV and a% = 0.0068 



In Fig. 1 we plot a typical solution, including also pb 
and /9 ra d- During radiation domination p v , p x -C Pradi 
then ip ~ const., so that DM exhibits the standard scal- 
ing, e -3 ^. Around matter-radiation equivalence the at- 
tractor starts to become effective and after some transient 
regime - possibly including a baryon-dominated epoch - 
it is reached at about the present epoch. In this scenario, 
the future will be characterized by values of £l v , and Q x 
closer and closer to the values in (5), with pb and p rac j 
diluting with the usual laws, a" 3 and a~ 4 , respectively. 

As will be discussed in [8] the growth of fluctuations in 
this scenario is compatible with the present data on CMB 
and large scale structure. In particular, it is remarkable 
that matter fluctuations grow even during accelerated ex- 
pansion, unlike in the usual A-dominated or quintessence 
scenario. 

2. VAMPs and DM abundance As we have discussed 
in the previous section, the late-time attractor forces 
the DM particle mass to change with time as to en- 
sure that DM and DE scale at the same, constant, ra- 
tio given by Eq. (5). In this section we show how this 
translates into an enhancement of the parameter space 
giving the observed energy densities. Alternative ideas 
to modify the standard freeze-out picture of CDM sup- 
pose a kineton-dominated expansion of the universe [9] or 
a low-reheating temperature after inflation [10]. Differ- 
ently from these approaches, where some new dynamics 
occurs at early epochs, in our proposal the standard ex- 
pansion history of the universe gets modified only at late 
times - after radiation-matter equality - when the field 
ip starts to roll. We assume that the DM particle x f aus 
into the category of Cold Dark Matter (CDM), that is 
X-particles cease annihilating during the evolution of the 
universe when they were non-relativistic ( in the case of 
relativistic DM, due to the fact that the DM abundance 
Y x is mass independent, the predictions of the VAMP 
and the standard scenarios are the same) . This amounts 
to assuming that Xf = M x (iff) /Tf > 1, where tpf is the 
value of the DE field at the freeze-out epoch and Tf is 
the freeze-out temperature. The present DM energy den- 
sity, Eq. (3), depends on the value of the physical x-mass 
today M£ h = M x (<p ) = Me~ A</5 ° (where (p is the value 
of the DE field today) and on the comoving number den- 
sity Y x . The latter depends, among other parameters, 
on the x mass at the earlier epoch of freeze-out, i.e. on 
M x (tp f ) = Me'^f = MP h e AA ^, where Aip = <p — ipf. 

In order to be definite, we will consider the case of sta- 
ble x-particles annihilating into light states by exchang- 
ing an heavier particle S of mass Ms, XX ~ * light states. 
To estimate xj we make use of the criterion r ann (a;/) ~ 
H(xf), where r ann = n^, q (<7 antl |i;|) and (<7 ann is the 
thermally-averaged annihilation cross-section times ve- 
locity of the DM particles. We also suppose that the DM 
particles annihilate through a p-wave suppression chan- 
nel. This happens - for instance ~ if x is a Majorana 
fcrmion. 
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Under these assumptions, the thermally-averaged an- 
nihilation cross section at freeze-out is given by 
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(8) 



where as is a combination of coupling constants 
measuring the strengths of the interactions. The 
standard freeze-out picture [7] gives the number 
abundance at late times Y X (T ) ~ Y x (Tf) = 

(7.Qx 2 f /gl /2 {T f )M p a a M x {ip f )), where T ~ 2.75 K is 
the present-day photon temperature and g* counts the 
relativistic degrees of freedom in the plasma at the freeze- 
out epoch. The energy density of the DM particles today 
reads 



£l x {to)h 2 = ce 



e 2\A v ( M ph) 
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(9) 



with c ~ 1.46 x lir w a' s 2 {x)lg 1 J 2 (T f ))GeSr- 2 . The 
standard ACDM scenario is recovered by setting Aip = 
or, equivalently, A = j3 = 0. In this case, V plays the role 
of the cosmological constant, M x = AfP h at any epoch, 
and constant £l x (to)h 2 contours in the (M x — Ms) plane 
are those corresponding to <f> = 1 in Fig. 2, where the 
region Ms < M x has been excluded. It is worth while 
stressing that, once the acceptable range for £l x (to)h 2 has 
been chosen, the DE abundance has to be fine tuned to 
give a flat universe (at, say, ten percent accuracy). This 
translates into a constraint on V which is nothing but 
the cosmic coincidence problem of the ACDM scenario. 
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FIG. 2. Regions with 0.1 < fl x h 2 < 0.3 and Ms > M x 
for fixed values of Aip, see Eq. (10). The grey region repre- 
sents the parameter space available in the standard scenario 
Aip — 0. The red+grey region represents the total parameter 
space available in the VAMP scenario. 



In the VAMP scenario the situation changes dramat- 
ically. Keeping the physical abundance, fl x (t )h 2 , fixed 
and assuming the attractor has been reached by today, 
the circular contour for Aip = is deformed to a ellipse 
of the same area and eccentricity given by e* Aip , as in 
Fig. 2. On the same contour the potential energy of the 
ip field is also constant, Ve" Vo , so we can express Aip in 
terms of V as 



1 , V 1 . _ 
A^--log^--log$ 



(10) 



where V is the particular value of V such that Aip = 0. 

Compared to the ACDM case, in the VAMP case we 
can vary V by many orders of magnitude, see Fig. 2. 
In principle, any value of V would give an acceptable 
ellipse. In practice, extreme values would prevent the 
field ip from reaching the attractor before the present 
epoch, thus deforming the corresponding contour from 
an elliptical shape and reducing its area. However, as 
we have checked numerically, the allowed values for V 
typically span many orders of magnitude, thus removing 
any fine tuning on this parameter. 

The most striking feature is that observationally al- 
lowed values of £l x (to)h 2 can be obtained for very large 
masses of the scalar particle S. In the traditional ACDM 
scenario large scalar masses would lead to a suppression 
of the annihilation cross-section and therefore to a over- 
abundance of DM particles today. On the other hand, in 
the case in which the x-mass depends on the DE field, 
the initially large value of the %-energy density at the 
freeze-out is adiabatically reduced not only by the ex- 
pansion of the universe, but also by the fact that, due to 
the running of the DE field, the x-mass is smaller today 
than at freeze-out. For consistency, we have to impose 
in this case that x has always been lighter than S, which 
is ensured by the constraint M x (ipj) < Ms, giving the 
straight dashed lines in Fig. 2. 

Even more remarkable from a phenomenological point 
of view is the drastic enhancement of the allowed param- 
eter space in the (M^ h — Ms) plane. Varying V, the 
allowed regions span an area delimited by an hyperbole, 
which contains values for Ms much larger than in the 
standard scenario. 

Points above the hyperbole give still the good DM/DE 
ratio, but fail in reproducing their absolute value, that 
is, the Hubble parameter today. 

4- An illustrative example: super symmetric dark mat- 
ter. Of the many CDM candidates, the best motivated 
from a particle physics point of view is the supersymmet- 
ric neutralino if the latter is the lightest supersymmetric 
particle (LSP) [11]. Most of the supersymmetric mod- 
els obtained from supergravity usually predict that the 
LSP is an almost pure f?-ino. Therefore, to give an il- 
lustrative example of our previous findings, we consider 
the case in which the LSP is a B-'mo and its mass Mb 
depends exponentially upon the time-dependent DE field 
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<j), M B = M B exp(-A^). 

Being the i?-ino the supersymmetric partner of the 
U(1)y gauge boson, it does not directly contribute at 
the one- loop level to the running of gauge couplings. This 
means that the time-dependence of the electromagnetic 
coupling is induced by the bino mass only through tiny 
quantities parametrizing how much the LSP composition 
depends on the Wino and Higgsinos. A full parameter 
space analysis regarding this specific point will be pre- 
sented in [12]. In the early universe, i?-inos mainly an- 
nihilate into fcrmion pairs through i-channel exchange of 
squarks and sleptons. Exceptions occur only for patho- 
logical situations in which there is a resonant s-channel 
exchange of Z° or a Higgs boson. Because of the large hy- 
percharge of the right-handed electron and the expected 
lightness of sleptons compared to squarks, it is often a 
good approximation to include in the annihilation cross 
section only the exchange of the right-handed sleptons. 

Summing over three slepton degenerate families with 
mass Mf R , the thermally-averaged £>-ino annihilation 
cross section times velocity has the same form as eq. (8), 
with Ms and M x replaced by Mp R and M B , respectively, 
and a| = 2Aira 2 / cos 4 W ~ 6.7 x 1(T 3 . 

In the traditional case in which the mass of the B-mo 
is not dependent upon the DE field, Mb = M B , cosmo- 
logical considerations give an upper bound to the i?-ino 
mass. Indeed, the requirement that charged particles are 
not the LSP implies Mt R > Mb- The minimum allowed 
B-mo relic abundance corresponds to the maximum anni- 
hilation cross section and therefore to the minimum M^ R . 
Setting Mi R = Mb in the expression for 0, one obtains 
an upper bound on the B-mo mass of about 300 GeV (for 
flh 2 < 0.3) [13]. This bound can be weakened in several 
ways. For example through the presence of resonant s- 
channel annihilations, once a small Higgsino admixture 
is introduced. Furthermore, as emphasized in Ref. [14], 
whenever the sleptons and the i?-ino become degener- 
ate in mass within about 10-20%, one cannot ignore the 
effects of coannihilation. These effects can modify signif- 
icantly the -B-ino relic abundance, because annihilation 
channels involving the charged sleptons have large cross 
sections which are not p-wave suppressed. Indeed, even 
in the case of the constrained model, the previous limit 
on the 5-ino mass can be relaxed to about 600 GeV 
[14]. Coannihilation effects do not significantly modify 
the bound on the slepton mass for a fixed value of Mb 
(as long as it is not too close to Mt R ). 

On the other hand, these bounds on the slepton and 
B-ino masses can be drastically modified if the mass of 
the B-mo depends upon the DE field. These effects are 
illustrated in Fig. 2 at different values of the parameter 
V . As expected from the findings of the previous sec- 
tion, observationally allowed values of f2.s(to) h 2 can be 
obtained for very large slepton masses for given values of 
V: the smaller V is, the larger the slepton mass can be. 
Furthermore, allowing V to vary enhances drastically the 
allowed parameter space in the (M B h — Mg R ) plane. We 



leave a more complete analysis of the parameter space of 
the supersymmetric models for a future publication [12]. 
It is intriguing that a coupling between the dark energy 
and the dark matter sectors may provide not only a nat- 
ural solution to the coincidence problem but also a dras- 
tic modification of the predicted abundance of CDM in 
terms of the particle physics parameters. 
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